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Abstract. Elastic phase transitions of crystals and phase transitions 
whose order parameter couples linearly to elastic degrees of freedom 
are reviewed with particular focus on instabilities at zero temperature. 
A characteristic feature of these transitions is the suppression of critical 
fluctuations by long-range shear forces. As a consequence, at an elas¬ 
tic crystal symmetry-breaking quantum phase transition the phonon 
velocity vanishes only along certain crystallographic directions giving 
rise to critical phonon thermodynamics described by a stable Gaussian 
fixed point. At an isostructural solid-solid quantum critical end point, 
on the other hand, the complete suppression of critical fluctuations re¬ 
sults in true mean-field critical behavior without a diverging correlation 
length. Whenever an order parameter couples bilinearly to the strain 
tensor, the critical properties are eventually governed by critical crystal 
elasticity. This is, for example, the case for quantum critical metamag¬ 
netism but also for the classical critical Mott end point at finite T. We 
discuss and compare the solid-solid end points expected close to the 
Mott transition in V 2 O 3 and k-(BEDT-TTF) 2 X. 


1 Introduction 


In most quantum or classical phase transitions in solids, the underlying lattice is 
affected by criticality. The change of the lattice constant measured in a thermal 
expansion experiment is, for example, often a very sensitive probe of criticality. In 
this review, we will focus on situations where either the lattice itself becomes critical or 
strongly affects criticality. The ratio of thermal expansion, a = {l/V){dV/dT)p, and 
specific heat, Cp = T{dS/dT)p, turns out to be a useful, unbiased tool to identify and 
classify zero temperature quantum phase transitions in general, that can be tuned 
with pressure. It can be identified with the relative change of temperature upon 
adiabatically varying the pressure 


r = 


1 

V T 



a 



( 1 ) 


with the molar volume Vm, and it is closely related to the parameter originally in¬ 
troduced by Griineisen [T]. The Griineisen ratio F is predicted to diverge with a 
characteristic power-law at a quantum critical point [2]. At the critical pressure Pc, it 
diverges as a function of temperature as T ~ with an exponent determined 
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by the product of the correlation length exponent v and the dynamical exponent 2 . 
Moreover, in the low-temperature limit the divergence as a function of pressure is 
universal, F = 


y ^ , in the sense that its prefactor G is solely determined by 
critical exponents. The Griineisen parameter thus not only provides a necessary crite¬ 
rion for the existence of a quantum critical point but also yields important information 
about its universality class by identifying certain combinations of critical exponents. 
Moreover, F also changes sign close to a quantum critical point, and the sign changes 
reveal the locations of entropy acummulation [3] . Measurements of F mm and its 
magnetic counterpart, the magnetocaloric effect |8I9I10I11I12I13) . have proven use¬ 
ful in the investigation of various quantum critical systems, see Refs. [14115116) for 
reviews. 

However, the theory of Ref. [3] implicitly assumes that the crystal lattice couples 
only perturbatively to the critical system and can be considered as a non-invasive 
probe that does not modify the critical properties itself. This assumption is often jus¬ 
tified. Nevertheless, there are exceptions where the crystal lattice and critical degrees 
of freedom are instead strongly coupled and the assumptions of Ref. [2] break down. 
Such quantum phase transitions that are governed by crystal elasticity are at the 
focus of this short review. We will argue that the non-perturbative coupling to the 
crystal lattice has a profound effect and fundamentally changes the critical behavior. 
The reason for this is that the rigidity of the crystalline lattice can induce long-ranged 
forces which suppress critical fluctuations, reminiscent of the suppression of Goldstone 
modes by long-ranged gauge interactions within the Higgs mechanism. Gonsider, for 
example, a large fluctuating domain of size L close to a critical point where the lattice 
constant within this domain is changed by a factor 1-l-e with e <?; 1. If the product Le 
is much larger than a lattice constant a, the crystal becomes strongly strained. As a 
result, the corresponding energy cost can suppress the formation of such domains that 
exceed a certain size. It can, however, happen that the whole crystal deforms without 
such large scale fluctuations in a transition described asymptotically by mean-field 
theory. 

For classical phase transitions, the coupling of critical fluctuations to the elastic 
degrees of freedom, following early work by Rice m and Domb [18] , were extensively 
studied in the 1970ies |19l20l21l22l2:il24l25l26l27l28l29l:-;()l:H| . Some of the results 


presented here are straightforward generalizations of these previous works to the case 
of quantum phase transitions. We start in section |1.1| with a short introduction to 
elasticity theory and recapitulate in section [l. 2 | the most important elastic couplings 
to the critical fluctuations. Sectionreviews Ref. [32] and discusses quantum critical 
elasticity focussing in particular on spontaneous crystal symmetry-breaking quantum 
phase transitions and isostructural solid-solid quantum critical end points (QGEPs). 
In section]^ we discuss quantum critical metamagnetism and argue that it is generally 
preempted by critical elasticity providing an example of a fluctuation-induced solid- 
solid QGEP. Finally, in sectionj^the critical properties of the Mott end point at finite 
temperature is reviewed, which is in fact an example of a classical solid-solid end 
point. We distinguish between a strong and a weak non-perturbative elastic coupling 
scenario that might be relevant for V 2 O 3 and k-(BEDT-TTF) 2 A, respectively. Some 
of the results for the Mott end point have been already presented in Ref. |33|. Part 
of this review is also based on Ref. [34] . 


1.1 Basics of elasticity theory 

The fundamental quantities of elasticity theory are introduced in the following, which 
are required for the subsequent discussion. For an in-depth introduction we refer the 
reader to standard textbooks [33I3H| . 
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The static elastic properties of a crystal are described by the following potential 
for the strain tensor field Sij 


Telast — (2) 

where Cijki is the elastic constant tensor. It is convenient to decompose the strain 
field into a macroscopic part Eij and a part that carries finite momentum, 

Eij (r) = ^ (diUj (r) + dju, (r)^. (3) 

The latter is related to the phonon displacement field Ui. The low-energy phonon 
Lagrangian reads 


PrEi \2 1 


■^phonon — 


(4) 


with the mass density p. The phonon energy-momentum dispersion follows from the 
Lagrangian Q and is given by puj'^ — qiCijkiqk = 0. In particular, the phonon velocities 
are determined by the eigenvalues of the dynamical matrix Mj;(q) = qtCijkiqk- On 
the other hand, the macroscopic strain is governed by the potential 

~ '^^ijC'ijklEkl (JijEij (5) 


where we also allowed for the presence of a macroscopic stress tensor aij. In case of 
an applied hydrostatic pressure, p, the stress tensor is diagonal aij = —pSij. For a 
uniaxial pressure, say, along the x-direction, aij = —px^ix^jx etc. The macroscopic 
strain also possesses kinetic energy, and the boundary conditions couple its equation 
of motion to the one of the phonon modes. Both effects do not contribute in the 
thermodynamic limit and will therefore be neglected in the following. 

Generally, the strain Eki can be decomposed into irreducible representations of 
the crystal class. Depending on the degeneracy of the corresponding eigenvalue of 
the elastic constant matrix, these representations can be either singlets, doublets or 
triplets. For example, a cubic crystal possesses a singlet Ej^x^ + Eyy + E^z, a doublet 
{Exx - Eyy, 2Ezz - Exx - Eyy) uud u triplet {Eyz, E^z, Exy) with eigenvalues of the 
elastic constant matrix in Voigt notation Cn -|-2Ci2, Cu — C 12 , and C 44 , respectively. 
Macroscopic stability of the crystal requires that each of these eigenvalues are positive, 
i.e., in general the tensor Cijki should be positive definite so that the potential ([^ is 
stable. 


1.2 Elastic coupling to critical fluctuations 

We now consider the coupling of these elastic degrees of freedom to the critical fluc¬ 
tuations of a second-order phase transition. In order to be concrete, let us assume 
that this transition is described by the following (/("^-potential in terms of a real singlet 
order parameter (f> 


Vcr = ^(t>^ + - h(j) (6) 

where u > 0 and the critical point obtains for r = 0 and h = 0. The order parameter 
(j) might represent some magnetic, electronic, or other degrees of freedom. The crystal 
lattice that houses a critical subsystem will unavoidable couple to the fluctuations of 
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the order parameter ij). Let us consider the most important couplings that are linear 
in the strain field. As we consider here only a singlet order parameter (j), this coupling 
can only involve a certain singlet of the crystal group that we denote by e, and the 
interaction can be written in the form 

Vint = (7) 

Whereas the coupling 72 is expected to be finite, the bilinear coupling 71 is often 
forbidden by symmetry. For example, for an Ising ferromagnet where (j) corresponds 
to the Ising magnetization the coupling 71 must vanish due to time-reversal symme¬ 
try. Similarly, order parameters that break translational symmetries cannot couple 
linearly to the strain. The coupling 71 is however finite, for example, at end points 
of lines of first-order metamagnetic or Mott transitions, or at an electronic nematic 
transition. 


1.2.1 Perturbative treatment of the elastic coupling 

In the perturbative regime, we can neglect the phonon contributions and consider 
only the macroscopic strain singlet, E, in Eq. 0 . For a cubic crystal, for example, 
the elastic part then reduces to 

Vei = ^l 2 E(j)'^ - 71 -b y - pE ( 8 ) 

where K is the corresponding modulus, i.e., K = (Cn + 2Ci2)/3, and p is the hydro¬ 
static pressure. Perturbatively minimizing for the strain singlet one obtains E « pjK 
and the effective critical potential assumes the same form as Eq. 0 

Veff = - hip)4' (9) 

but with pressure-dependent parameters, r{p) = r + ^ipjK and /i(p) = h + ^ipjK. 
The tuning parameters, i.e., the prefactors of the most relevant operators in the 
critical theory, h and r, thus obtain an effective pressure dependence. 

Let us first discuss the case where 71 is forbidden and, furthermore, h = 0. The 
contribution to the free energy attributed to the quantum critical degrees of freedom 
can then be written in the scaling form 



r{p) 

T^Kvz) 


)■ 


( 10 ) 


where the scaling function / is assumed to possess the following asymptotic behavior 
/(cc) —>■ const, for X —?► 0 and /(x) ~ for |x| —?> 00 . This scaling form for the 

quantum critical free energy with a pressure dependent tuning parameter, r, was used 
in Ref. 0^for the discussion of the Griineisen parameter. However, the critical free 
energy ( |l0| ) also predicts a critical renormalization of the bulk modulus. At T = 0 it 
reduces to J^cr oc and the renormalization of the bulk modulus reads 


( 11 ) 


This renormalization will be perturbative as long as 

v{d + z) — 2 > 0. 


( 12 ) 
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This is just the generalization of the corresponding result for classical criticality where 
the elastic coupling remains perturbative if the specific heat exponent is negative 
|m 8 l‘^l| . i.e., a = 2 — i/d < 0, from which the criterion for quantum criticality 
( |12[ ) obtains after replacing the dimension d with the enhanced dimension d + z. This 
criterion simply arises from the requirement that the perturbative renormalizations of 
the elastic constants due to the correlator (</>^^^) of the order parameter remains non¬ 
singular. Although there might be exceptions for effectively low-dimensional systems, 
see e.g. Ref. the criterion ( |12[ ) for quantum criticality is usually fulfilled. In fact, 
the theory of Ref. [2] for the Griineisen parameter implicitly assumes that the criterion 
([T^ is obeyed. 


1.2.2 Non-perturbative elastic coupling to critical degrees of freedom 


The elastic coupling 72 becomes non-perturbative if the criterion (12) is violated. On 


the other hand, the bilinear coupling 71 will always be non-perturbative. In the latter 
case, the corresponding perturbative correction to the bulk modulus is attributed to 
the correlator {(pcj)) of the order parameter that is singular by definition. Note that this 
also applies to situations where the order parameter (j) itself is a bilinear in fermionic 
fields and the critical theory must be formulated in terms of fermionic degrees of 
freedom. Whenever a coupling 71 is present it will lead to a strong entanglement be¬ 
tween the critical and elastic degrees of freedom so that a perturbative treatment of 
the elastic coupling unavoidably breaks down. If no first-order transition is induced, 
one obtains instead an elastic critical point. The properties of such continuous elas¬ 
tic phase transitions is at the focus of the following chapters. In chapter we start 
with a general analysis of continuous elastic quantum phase transitions irrespective of 
their origin, i.e., whether they are genuine or are induced by other, critical degrees of 
freedom. The effective critical elastic theory is analysed, the quantum critical thermo¬ 
dynamics, and, in particular, the behavior of the Griineisen parameter is discussed. 
In chapter the particular example of metamagnetic quantum criticality is investi¬ 
gated that generically allows for a bilinear coupling 71 and is thus non-perturbative. 
Finally, in chapter]^ a corresponding example of classical criticality is also discussed 
where we review the issue of the universality class of the critical end point at finite 
temperature for the Mott transition. 


2 Quantum critical elasticity 

The instabilities of a crystal lattice can be studied from the point of view of an 
effective critical theory in terms of the strain order parameter. Such an approach 
was chosen in Ref. [32] where the quantum critical signatures of continuous elastic 
quantum phase transitions were discussed building on previous work by Cowley [55] 
and Schwabl and collaborators |2Vl28l2Hl3()| . At such an instability an eigenvalue of 
the elastic constant matrix Cijki vanishes. The strain order parameter is identified as 
the corresponding eigenvector, and depending on the degeneracy of this eigenvalue it 
can either be a singlet, doublet or a triplet of the irreducible representations of the 
crystal class. 

Importantly, if such an eigenvalue of Cijki vanishes the phonon velocities in general 
remain finite. Mathematicaly this arises because the eigenvalues of Cijki are deter¬ 
mined from the 6 x 6 matrix of the elastic constant matrix in Voigt notation, Cp\, 
whereas the phonon velocities derive from the eigenvalues of the dynamical matrix 
= diCijkidk (see above). Physically, it reflects the fact that general deforma¬ 
tions of the crystal cannot be represented in terms of an acoustic phonon mode for 
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elastic transition 

constant 

strain 

type 

orthorhombic —>■ monoclinic 

C44 

£23 

I 

orthorhombic —>■ monoclinic 

C55 

£13 

I 

orthorhombic —>■ monoclinic 

C66 

£12 

I 

tetragonal —>■ orthorhombic 

Cll — C12 

£11 — £22 

I 

tetragonal —>■ orthorhombic 

C66 

£12 

I 

tetragonal —>■ mono- or triclinic 

C44 

(£23,£13) 

I-fll 

hexagonal —> mono- or triclinic 

C44 

(£23, £13) 

I-fll 


Table 1. Table from Ref. |32] based on results of Refs. | 26I27I29| listing possible second- 
order elastic phase transitions. The columns from left to right specify the crystal symmetry 
that is broken at the elastic transition, the eigenvalue of the elastic constant matrix in Voigt 
notation that vanishes at the transition, the strain order parameter and the type of the 
transition in the classification of Cowley [^. Certain modifications apply for tetragonal 
crystals with a finite component cie. 


fixed direction, q = q/|q|. Only for specific directions of the phonon wavevector q, 
the polarization of the phonons might induce the particular type of strain that is 
akin to the one of the strain order parameter so that the phonon velocity only van¬ 
ishes for this particular direction q in momentum space. As discussed by Cowley [13] 
and Folk et al. [27], there are elastic transitions where the phonon velocities either 
remain finite for all q, vanish for q pointing along certain directions or lying within 
specific two-dimensional planes of momentum space. The phonons thus become only 
critical in a m-dimensional subspace with m = 0,1, 2 corresponding to type 0, I or 
II, respectively, in the classification of Cowley [53]. As result of this restricted criti¬ 
cal subspace, even the classical critical theory is either at or above its upper critical 
dimension for m = 2 and m = 0,1, respectively [27l8T| . We will see below that the 
corresponding quantum phase transition is always above its upper critical dimension. 
In the following, we distinguish between elastic transitions that do and do not break 
a crystal symmetry. 


2.1 Spontaneous symmetry-breaking elastic transitions 

For elastic transitions that break a crystal symmetry the expectation value of the 
macroscopic strain order parameter E is zero in the symmetric undistorted phase 
and assumes a finite value in the symmetry-broken phase. Such a transition can be 
discontinuous or continuous, i.e., first or second order, respectively, depending on 
the properties of the Landau potential V{E) for the order parameter. Assuming that 
V{E) is analytic so that a Taylor expansion around E = 0 can be performed, the 
transition can be of second order whenever the cubic term in this expansion vanishes. 
According to Refs. [26127129] . this is the case for the transitions listed in Table 
With the exceptions listed in the last two rows in Table [^ most of these transitions 
are characterized by a singlet strain order parameter A, for which the potential reads 

V{E)='^E^ + ^E^ + aE, (13) 

where cr is the appropriate singlet of an applied stress, which explicitly breaks the 
symmetry E —>■ —A, and a cubic term E^ is absent. If the quartic coupling is positive 
M > 0, a second-order elastic quantum phase transition occurs at T = 0 and <7 = 0 
when the tuning parameter goes to zero r —^ 0"*". This is the case when the elastic 
constant listed in the second column of Table [^ vanishes. The tuning parameter 
r = r{T,p) in general depends on temperature, T, and hydrostatic pressure, p. The 
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I 

(i) 

(ii) 

Ccr- 

Ocr ~ 
OLcr / CcT 

/ 2 

rj^Z/2 

~ 1/T 

y3j,-l/2 

rp3^-3/2 

_ 1 1 
_ 6 P-Pc 


II 

(i) 

(ii) 

CXcr 

O^cT / Ccr 

T 

~ 1/T 

T^r-^ 

_ 1 1 

3 P — Pr 


Fig. 1. Left: Phase diagram for a symmetry-breaking elastic quantum phase transition. The 
parameter r(T = 0) can be tuned by pressure and vanishes for tuning the corresponding 
elastic constant to zero, see second column of Table Its temperature dependence r{T) 
determines the shape of the finite-temperature phase boundary. The thermodynamics at¬ 
tributed to the critical phonons shows a crossover at T ~ r resulting in two regimes (i) and 
(ii). Right: Critical phonon thermodynamics: specific heat Ccr, thermal expansion Ocr and 
Griineisen ratio in the regimes (i) and (ii) for pressure tuning r{T = 0) oc p — Pc', for type I 
(upper table) and II (lower table), respectively. 


latter effectively arises from anharmonicities in the elastic potential that couples 
the singlet E to the trace of the strain tensor, Vanha tr{eij}i?^. Minimizing the 
trace in the presence of a hydrostatic pressure results in a pressure dependence of 
K similarly as discussed in section EM The temperature dependence, on the other 
hand, is induced either by anharmonicities, i.e., phonon excitations or coupling to 
other degrees of freedom, e.g., electrons. 

The symmetry-breaking elastic transition are either of type I or II in the Cowley 
classification. The wavevector of the phonon can be decomposed, q = (p,k), into an 
m-dimensional soft component p and a (3 — m)-dimensional stiff component k with 
m = 1,2 for type I and II, respectively. Close to the elastic phase transition, the 
phonon dispersion then possesses the anisotropic form [28] 

^ rp^ -I- ap"^ -b 6k^ -b ... (14) 


with finite constants a and b, and the dots represent further terms that are less 
relevant. For the identification of scaling exponents, it is convenient to perform the 
substitution k^ —>■ k'^. This effectively corresponds to the introduction of an enha nced 
spatial dimensionality des = to -b 2{d — to) = 2d — m with d = 3. The dispersion (14) 
is then characterized by the scaling r ^ p^ and ~ p'^jk'^, which determines the 
correlation length exponent v = \/2 and the dynamical exponent z = 2, respectively. 


For classical elastic transitions, the upper critical dimension of the theory is 
obtained from the criterion deft = 4 which amounts to dpprp = 2 -b to/ 2 |27l31j . 
While classical transitions with to = 0,1 are above their upper critical dimension 
d > d)tprp for d = 3 allowing for a mean-field description, classical type II transitions 
with TO = 2 are at their upper critical dimension giving rise to logarithmic correc¬ 
tions to mean-field exponents. For elastic quantum phase transitions, on the other 
hand, the upper critical dimension is further reduced by the dynamical exponent, 
dqpT = 2 + to/2 — 0 = to/ 2. All elastic second-order quantum phase transition are 
therefore above their upper critical dimension for d = 3 and the phonon excitations 
can be treated perturbatively. 
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In particular, the contribution of the critical phonon excitations to the free energy 
can be written in the scaling form [5] 


J^cr=T 






(15) 


with vz = I and = 4 — m/2. An explicit calculation shows that the function 

/ possesses the asymptotics f{x) = const, for cc —>■ 0 and f{x) ^ _ ^-m /2 

for a; —> oo. From Eq. (15) the critical phonon thermodynamics is easily derived. At 
criticality r = 0, i.e., in regime (z) of Fig. the critical contribution to the phonon 
specific heat is given by Ccr ~ j ^ ^ 7^2 |.ypg j 

type II transitions, respectively. The volume thermal expansion, a, derives from the 
pressure dependence of the tuning parameter r so that Ocr y 2 -m /2 gt 7 - = Q. For 
the critical Griineisen ratio defined as //r = Ocr/C'cr follows 




1 


T^l{vz) 


(16) 


with nz = 1 in agreement with scaling considerations [2]- In the low-temperature 
regime (zz) of Fig. the divergence is universal 




m 1 
6 Vm{p-Pc) 


(17) 


with the prefactor m/ 6 , i.e., 1/6 and 1/3 for type I and II, respectively, where we used 
r(T = 0,p) oc p—pc with the critical pressure Pc- The critical phonon thermodynamics 
is summarized in the Tables of Fig. 

It should be noted, however, that the critical phonon thermodynamics vanishes 
with a relatively high power of temperature for T —>■ 0. As a consequence, it might be 
actually subleading compared to other non-critical contributions, e.g., due to gapless 
electronic degrees of freedom in metals. In such a case, the identification of the critical 
phonon parts, Ca and acr, and thus the analysis of the critical Griineisen ratio //.j. 
might be challenging. 


2.2 Isostructural elastic transitions 

Elastic transitions not listed in Table are generically first-order and exceptions 
require additional fine-tuning. For example, the Taylor expansion for the potential 
of the strain order parameter might allow for a cubic term whose prefactor can, 
however, be tuned to zero. Another example, on which we focus in the following, are 
isostructural transitions. The corresponding macroscopic strain order parameter, E, 
is a singlet representation that is invariant under all crystal symmetry operations. Its 
Landau potential generally includes all powers of E; after shifting the order parameter 
by a constant E —>■ E+Eq one can eliminate the cubic term and the potential assumes 
the form 


V{E) = ^E^ + ^E^ + aE, (18) 

with K 4 > 0. For AT < 0 a first-order phase transition arises where the expectation 
value of E jumps at ct = 0 as function of a. This corresponds to a transition between 
two different solids that possess the same crystal symmetries while they are distin¬ 
guished by different values of E. A second-order solid-solid quantum critical end point 
(QCEP) obtains when both parameters, K and cr, are tuned to zero at T = 0. Here, 
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Fig. 2. Phase diagram with a solid-solid quantum critical end point (QCEP). A line of first- 
order solid-solid quantum phase transitions meets a line of finite temperature second-order 
transitions at the QCEP and confine a surface (shaded area) of first-order transitions at 
finite T. The dashed coordinate system is defined by the parameters K and a, see Eq. ( |18[ ), 
that vanish at the QCEP. 


the expectation value of E changes in a critical manner, for example, as a function of 
a. The tuning might be achieved with the help of pressure p and an external field F, 
on which the parameters depend, K = K{F,p) and tr = a{F,p); the position of the 
QCEP within the {F,p) phase diagram is then identified by K{Fc,Pc) = a{Fc,Pc) = 0, 
see Fig. 


Minimizing the potential (181 one obtains 


E : 


a/K 

-\a/K4\^/hgnia) 


for \a\ < 
for |(t| » 


(19) 


In the latter limit at F = Fc, the strain depends on the pressure E ~ {p — PcY^^ in a 
non-linear fashion with mean-field exponent 5 = 3. This corresponds to a breakdown 
of Hooke’s law, which is a hallmark of solid-solid end points. 

Importantly, the solid-solid QCEP is characterized by the absence of critical mi¬ 
croscopic fluctuations as the velocity of sound remains finite for all directions in this 
case. In fact, all isostructural transitions are of type 0 in the Cowley classification 
and the phonon sector remains non-critical because of the high symmetry of the order 
parameter E. Consequently, the quantum phase transition associated with the solid- 
solid QCEP is a true mean-field transition and exactly described by Landau theory. 
In particular, there does not exist a diverging correlation length and the usual scaling 
hypothesis for the critical free energy does not apply. The critical thermodynamics is 
in fact non-universal and depends on other, non-critical degrees of freedom that in¬ 
duce an effective temperature dependence of the tuning parameters, K = K {F, p, T) 
and a = a{F,p, T). For example, assuming that a{F,p, T) = ao{F,p)+aT^, with, e.g., 
a; = 2 in a metal, it follows from minimizing the potential (18) a critical contribution 
to the specific heat Ccr and thermal expansion Ocr ~ 

for r —>■ 0 at iF = 0. In the same limit one obtains for the critical Griineisen ratio 


r„ 


1 


1 


T=0,F=Fa 


3(a; - 1) Vmip-Pc) 


( 20 ) 


that diverges with the inverse of ao{Fc,p) (x p — Pc in agreement with expectations 
but with a prefactor depending on the exponent x. 

As the solid-solid QCEP separates two solids with different volumes but with the 
same crystal symmetry, the discussion given above is in fact relevant for all quantum 
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critical end points in solids that are not associated with the breaking of a true sym¬ 
metry. Such critical end points are described with the help of a scalar order parameter 
(f> and are in the Ising universality class. However, in a crystal this order parameter 
will in general couple linearly to the strain singlet, i.e., Vint = so that the di¬ 

verging susceptibility {4>4>) will cause a isostructural instability with the concomitant 
change in the universality class. As will be shown below, in this case the instability 
of the macroscopic strain E preempts the critical endpoint of the (j) field in such a 
way that critical fluctuations are absent at the transition. This applies, for example, 
to the Kondo volume collapse transition at zero temperature, where the asymptotic 
critical behavior will be governed by a solid-solid QCEP |38l39j . It is also the case 
for the metamagnetic QCEP [40141) . where F corresponds to a magnetic field, which 
will be discussed in some detail in the next chapter. Similar reasoning applies to a 
metaelectric QCEP, where the Ising order parameter is identified with the variation 
in the longitudinal electric polarization, (j) = P — Per, which increases in a critical 
manner at some finite electric field. For the corresponding classical metaelectric end 
point, that exists at finite temperatures in KH 2 PO 4 , the absence of a diverging cor¬ 
relation length has been already demonstrated in Ref. |32]- Another famous classical 
analogue is the 7 — a transition of Ce [43144145] . In chapter we will discuss the 
Mott end point at finite T that is also predicted to be governed by critical elasticity 


m- 


3 Metamagnetic quantum critical end point 

Metamagnetism is associated with a superlinear rise of the magnetization M{H) at 
some finite magnetic field Hm [ISj. As a function of decreasing temperature such a 
metamagnetic crossover might develop into a first-order transition where the mag¬ 
netization instead jumps as a function of field. The resulting line of metamagnetic 
first-order transitions terminates in a critical end point {Tep,Hm) located at finite 
temperature Tep- If this end point temperature Pep can be suppressed by some addi¬ 
tional tuning, a metamagnetic QCEP arises for Tep —>■ 0. 

Metamagnetic quantum criticality was introduced in the context of Sr 3 Ru 207 in 
order to account for the anomalous behavior close to its metamagnetic field [47140149] . 
The universal critical thermodynamics close to such a metamagnetic QCEP has been 
discussed by Weickert et al. in Ref. [H] however without considering the feedback of 
the lattice to the critical behavior. Moreover, Zacharias and Garst [41] analyzed in 
detail the critical thermodynamics close to such a QCEP within the spin-fluctuations 
theory as proposed in Ref. m- In particular, it was pointed out that the metam¬ 
agnetic QCEP is intrinsically unstable with respect to a magnetoelastic coupling. 
The metamagnetic Ising order parameter is given by the change of the longitudinal 
magnetization close to the critical field (j) = M — {M )Ks Ising symmetry is 
only emergent and, as a consequence, a linear coupling 71 to the strain is generi- 
cally allowed. This will have the consequence that the metamagnetic QCEP will be 
preempted by an isostructural quantum phase transition. 

3.1 Solid-solid quantum critical end point induced by metamagnetic fluctuations 

The effective bosonic theory for quantum critical metamagnetism in a metal reads 


m 



(j) [r - (j) + ^(f)^ -h(l) + Cdy, 


■n 


( 21 ) 
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where u > 0^ h (x H — Hm measures the distance to the critical field, and the 
metamagnetic quantum critical end point obtains for r = 0. The dynamical part 


‘5dyn — 


dr I dr £dyn = 7^ XI ^ 


CJrik 


( 22 ) 


is governed by the Landau-damping, i.e., damping of the order parameter by particle- 
hole pairs in the metal. The resulting thermodynamics was discussed in detail in 


Refs, mm- 

A linear coupling 71 of the order parameter to the strain field, see Eq. 0 , is 
generally allowed and will modify the quantum critical properties. An exhaustive dis¬ 
cussion of the theory (211 linearly coupled to the crystal lattice including an analysis 
of phonon degrees of freedom can be found in Ref. [31]. As the resulting isostructural 
transition will be of type 0 where the phonons remain non-critical, see section | 2 . 2 [ 
we can limit ourselves here to a coupling of the order parameter to the macroscopic 
singlet, E, of the strain tensor only. Effectively, the coupling 71 can then be absorbed 
into a shift of the tuning parameter h —> h + ^iE. Integrating out the order parameter 
fluctuations one then obtains the effective potential for the macroscopic strain E, 


Ves{E) = ^E^ + J-„(r, h + 7 iA, T) 


(23) 


where Kq is the corresponding modulus and Tcr(r,h,T) is the critical free energy of 
Eq. (21 1 . The magnetoelastic coupling gives rise to a renormalization of the modulus 
of the form 


K = alVeff(A) ^Ko- jfxcr (24) 

where Xcr = (</></') = —5^-Acr is the singular part of the differential susceptibility. 
At h = 0 and zero temperature this susceptibility is simply given by Xcr = 1/?' 
for r > 0. Consequently, before the metamagnetic QCEP is reached for r = 0 the 
renormalized bulk modulus vanishes at r* = > 0. The metamagnetic critical 

point is effectively preempted by an isostructural QCEP that is realized for r* = 

lUKo- 

Using the results for Ecr{r, h,T) of Ref. [41], one finds for the effective potential 
at r = r* the explicit form close to the isostructural QCEP 


Vqcep(A) = + Ip A" + aoE 


{h -I- 71A)^ 
2R* 


fo{r*,T) 


(25) 


where we introduced the corresponding singlet stress Uo that couples to E and also 
explicitly added a quartic term with coefficient K 4 that stabilizes the potential (that is 
also renormalized by SK 4 ~ ujf/R*^). The coefficient R* parameterizes the effective 
stiffness, Xcr = 1 /R*, of the metamagnetic potential at r = r*. 


R* =r* 


riT^, 

r 2 


i?* < £2/3 
R* > £2/3 


(26) 


The remaining free energy independent of E is given by 


/o(r*,£) = 


r*<£2/3 

-f6 T2r*'^, 


(27) 
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Fig. 3. Density plot of the metamagnetic correlation function x(q) of Eh- (30 1 in the limit 
|q| —>■ 0 within the {qx, qy) plane (a) in the absence, 71 = 0 and (b) in the presence, 71 ^ 0 , of 
a magnetoelastic coupling for a certain choice of parameters (r = 2,'y\IKQ = 1 . 3 , 7 i/(Cii — 
C 12 ) ~ 0.7, and 71 /C 44 = 0.5). For 71 0 the metamagnetic fluctuations inherit the spatial 

anisotropy of the phonon spectrum with a characteristic four-fold symmetry for the cubic 
crystal lattice. 


The coefficients ft and are explicitly listed in Ref. [41]. The potential (25) has 
the same form as the one given in Eq. (18) with parameters K = Kq — 'yf /Tr and 
cr = (To — "fih/ R*. 

A hallmark of the isostructural QCEP is the non-linear response of the crystal 
lattice with respect to an externally applied stress ctq corresponding to the breakdown 
of Hooke’s law. Minimization of the potential at T = 0 and h = 0 yields E cx CTq' . On 
the other hand, the metamagnetic fluctuations that drive the isostructural transition 
are reflected in a crossover at a characteristic temperature T ^ The specific 

heat for example, that follows from /cr(r*,T) in Eq. ( |2^ , shows an anomalous T- 
dependence, C ~ Ts at high temperatures, T 3> whereas it is of Fermi-liquid 

type, C ~ Tr*^~, for T Apart from the contribution to thermodynamics 

deriving from the underlying metamagnetic fluctuations, /o(r*,r), there will be also 
additional T-dependences of the bare parameters. In particular, the bare stress aq 
will possess an intrinsic temperature dependence in a metal, <Jq{T) = (To(0) -faT^, 
resulting in a thermal expansion of Fermi-liquid type, dxE = 8 x^ 0 /Kq ^ T, in 
the absence of criticality. Close to the solid-solid QCEP, however, this will result 
in an additional, critical contribution to the specific heat. For example, at h = 0 
and at lowest temperatures, Ccr ~ Td'^\aq{T)\^/^ ~ |(To( 0)|^/^T. Similarly, for the 
thermal expansion we obtain in the same limit Ocr |o’o(0)|~^^^T reproducing the 
result for the critical Griineisen parameter of Eq. ( |20| ) with x = 2. Note, however, 
that the critical contributions are here subleading compared to the ones deriving 
from fq{r*,T). In order to obtain the expected scaling for the critical Griineisen 
parameter /kr = ctcr/Ccr a careful subtraction of the leading non-critical contributions 
are necessary. This limits the practicability of the Griineisen analysis in the present 
case. 


3.2 Deviation of metamagnetic correlations from Ornstein-Zernicke form 


Whereas the phonons remain non-critical at the isostructural QGEP they nevertheless 
renormalize the metamagnetic fluctuation spectrum. In order to simplify the following 
discussion, we again assume a cubic crystal where the metamagnetic fluctuations 
couple linearly to the trace of the strain tensor. The quadratic part of the theory (21) 
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in the presence of an elastic coupling to the phonon modes Ui, see Eq. ([^, reads 

[r - V^] (t> - 'yi(j)d,Ui + ^d.UjCijkidkUi- (28) 

The static correlator of the metamagnetic fluctuations x(q) = for 7 ^ = 0 

possesses a standard Ornstein-Zernicke form 

Xo(q) = —^ (29) 

r + 

at /i = 0 and in lowest order in u. A finite linear magnetoelastic coupling 71 leads to 
a modification of the momentum dependence, 


x(q) 


xo(q) 

1 - 7?Xo(q)%(Af-i)ijqj 


(30) 


where M ^ is the inverse of the dynamical matrix Mji(q) = qtCijkiqk, see section 
0 For a cubic lattice the dynamical matrix reads explicitly 

Mij = {C 12 + Ci4)qiqj + + (Cn — C 12 — 2C4^4)qf^5ij. (31) 

where the last term should here not be summed over repeated indices. The bare 
modulus of the previous section is here identified with Kq = |(C'ii + 2Ci2). Due to 
the phonon renormalization, the correlation function x(q) depends on the orientation 
of the momentum q even in the limit |q| —>■ 0. This is illustrated in Fig.j^ The inverse 
l/x(q) assumes for |q| —> 0 the limiting values 

[ ^ - K,+ 2 {ct-C,,)/s for gl I (100) 
limoX-i(q) = \r- ^11 (HO) (32) 

i ^ - Ko+lc44 forgll(Hl). 

This explicit calculation also nicely demonstrates that the fluctuations indeed remain 
gapped at the isostructural QCEP that is realized at r = r* = yJ/ATo. Due to the 
stability conditions Cn — C 12 > 0 and Cn > 0, the gap lim|q|_,,o x~^(q) remains 
finite for all directions of q even at the isostructural QCEP. 

Such an asymmetry of the magnetic fluctuations is generically expected close to a 
metamagnetic end point and could, in principle, be detected with the help of neutron 
scattering. A similar asymmetry has been discussed and experimentally verified in 
the case of critical piezoelectric ferroelectrics [50] . 


4 Critical Mott end point at finite temperature 

Strongly correlated materials with a conduction band close to half-filling might exhibit 
a Mott metal-insulator transition as a function of doping or applied pressure m- 
When all lattice sites are on average occupied by a single electron, their motion is 
suppressed by a strong Coulomb repulsion that energetically disfavours to have a site 
doubly occupied giving rise to a Mott insulating state. The application of pressure, 
however, enhances the overlap integral of electron wavefunctions on adjacent lattice 
sites and thus the kinetic energy. For sufficiently large pressures the kinetic energy 
prevails resulting in a Mott metal-insulator transition at some critical pressure Pc- 
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Pc pressure p 


Fig. 4. Phase diagram of a system with a line of first-order Mott metal-insulator tran¬ 
sitions that terminates in a second-order critical end point at (pc,Tc)- Close to this end 
point elastic criticality prevails indicated by the yellow shaded regime. The criticality of the 
Mott end point therefore belongs to the universality class of a classical solid-solid end point 
characterized by mean-field Landau criticality |33| . 


This transition is first-order at low temperatures resulting in a line of first-order 
transitions in the (T, p) temperature-pressure plane. This line ends at a second-order 
critical end point at a finite T^. For temperatures T > there remains only a 
crossover so that the system can be smoothly transformed from an insulator to a 
metal by varying pressure and temperature, see Fig. The topology of the phase di¬ 
agram resembles the one of the liquid-gas transition and, correspondingly, it has been 
discussed theoretically and experimentally whether the Mott end point also belongs 
to the Ising universality class |52l53l54l55l5fil57l58l59lfi0l61162163184165166167168) . 

However, we will argue that this is in fact not the case. As the Mott transition 
is susceptible to pressure tuning its conjugate, i.e., the volume will exhibit a jump 
across the line of first-order transitions. Correspondingly, the black line in the phase 
diagram of Fig. |^does not only separate the Mott insulator from the metal but also 
two isostructural solids with different volumes. Consequently, its termination point is 
not in the liquid-gas, i.e., Ising universality but it is rather a classical critical solid- 
solid end point characterized by mean-field Landau criticality as already pointed out 
in early work by Jayaraman et al. |54j . 

In this section, we will discuss the critical behavior close to the finite-temperature 
Mott end point in some detail, especially as in this case quantitative estimates of the 
importance of the lattice coupling are available. We will not consider the correspond¬ 
ing Mott quantum phase transition at T = 0 here but, in contrast to the previous 
sections, strictly limit ourselves to the classical transition at finite Tc > 0. We will 
take it for granted here that the finite T-end point of the line of first-order Mott 
metal-insulator transitions on an incompressible lattice is indeed described by the 
classical Ising model. 


(33) 

with the Ising order parameter (p and pressure and temperature dependent tuning 
parameters h = h{j),T) and r = r{p,T). The end point of the incompressible lattice 
obtains for pressure and temperature values for which h = r = 0. The interaction 
u gives rise to a Ginzburg scale [5S] that defines a crossover between Mott mean- 
field behavior away from the end point, that is described by Landau theory, and 
non-trivial critical behavior ascribed to the Ising universality class. The coupling 
of (j) to the crystal lattice, Eqs. Q, eventually gives rise to another crossover to a 
regime governed by critical elasticity. Depending on the relative strength of the elastic 
coupling, one can distinguish two scenarios that we will discuss in the following. 






Will be inserted by the editor 


15 




Fig. 5. Sketch of two scenarios for crossovers upon approaching the Mott end point, (a) For 
relatively strong elastic coupling, there is a single crossover into a regime (yellow shaded) 
governed by elastic criticality, (b) For a relatively weak elastic coupling, first a crossover 
occurs into a regime (blue shaded) governed by the non-trivial Ising universality class, and 
afterwards a second crossover takes place into a regime (yellow shaded) controlled by elas¬ 
tic criticality. While scenario (a) applies to V 2 O 3 , scenario (b) is applicable to k-(BEDT- 
TTF)2X. 


4.1 Strong elastic coupling: Mott-Landau-to-elastic-Landau crossover 


For a relatively strong elastic coupling, there is just a single crossover from the mean- 
field Mott-Landau regime of the theory (33) to the regime in the immediate vicinity 
of the solid-solid end point that is governed by crystal elasticity, see Fig. [^a). This 
scenario can be captured within simple mean-field theory and is described by the 
potential. 


v(</), E) = -H+ - ^,Ec4 + ^E^ -{p- p,)E. (34) 


The curvature tensor dadpV{4>, E) with a, (3 = E, cj) at h = 0 and p = p^ toi (j) = 0 
possesses a vanishing eigenvalue at r* = ^\IKq. If the end point is approached for 
p = Pc and /i = 0 as a function of temperature, the Mott end point located at r = 0 
for the incompressible lattice is preempted by an isostructural solid-solid end point at 
r = r* > 0. For the strong elastic coupling scenario, the scale r* in addition preempts 
the Ginzburg crossover present for the Mott transition on the incompressible lattice. 
As a result, the critical properties are always captured by Landau mean-held theory. 

Nevertheless, one could dehne a crossover between a Mott-Landau regime where 
the elastic coupling is still perturbative and a elastic Landau regime where Hooke’s 
law breaks down. Solving the classical equations, d^V{(j), E) = 0 and dEVi4>: E) = 0, 
for the macroscopic strain, say, at criticality r = r*, h = 0, setting 72 = 0 for 
simplicity, one hnds that for pressure values \p — Pc\ S> ^p* with Ap* = "fijy/Kou, 
the response is linear, E « {p — pc)/Ko but becomes non-linear, E ^ \p — 
for Ip — Pel ^ Ap*. The pressure scale Ap* dehnes the width of the Landau critical 
regime along the pressure axis that is governed by crystal elasticity (orange shaded 
regime in Fig. [^a)). 

The scenario of a strong elastic coupling with a single crossover between a critical 
Mott-Landau to a elastic Landau mean-held regime might be applicable to V 2 O 3 . 
In Ref. [SB] the critical behavior of Cr-doped V 2 O 3 , i.e., (Vo.gsgCi'o. 011 ) 03 , close to 
its end point was experimentally investigated with the help of conductivity measure¬ 
ments. Its behavior could be explained with the help of mean-held exponents in a 
large range of pressures and temperatures except in the immediate vicinity of the 
critical end point. It has been speculated [SB] that the deviations observed very close 
to the critical pressure in a narrow range Ap/pc ~ 3%o indicate a Ginzburg crossover 
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Fig. 6 . (a) The second derivative —d'^T in dimensionless units at the critical pressure, Pc- 
It shows a pure mean-field jump (solid line) as a function of r preempting the logarithmic 
Ising singularity (dashed line), (b) Macroscopic strain as a function oi p — pc for various 
temperatures T, see text. 


to non-trivial Ising exponents. However, detailed measurements of various sound ve¬ 
locities by Nichols et al. [70] have revealed that the undoped compound V 2 O 3 , that 
translates to an effective distance from the critical end point of Apjp^ « 1.7% in 
pressure units, possesses already a substantial softening of the corresponding bulk 
modulus K. This modulus K denoted as A 4 in Ref. [70] has a deep minimum around 
540 K where it drops to around 26% of its value at room temperature. From our 
considerations above follows the expectation that the value of this A 4 modulus indeed 
vanishes exactly at the critical end point. The fact that it is so strongly reduced in 
the pure compound indicates that V 2 O 3 already accesses the elastic critical regime 
(yellow shaded in Fig. [^a)) as a function of temperature. In turn, it follows that 
most of the experiment of Ref. [56j where the conductivity at the critical tempera¬ 
ture was tuned with the help of pressure is located in this elastic critical regime. This 
suggests that the end point in this material is always governed by Landau mean-field 
theory and, in particular, that the deviations mentioned above are not caused by a 
Ginzburg crossover but are rather attributable to some other origin as, for example, 
the unc ert ainty in determining the precise form of scaling variables like r,h 01 p — Pc 
in Eq. (ISdl). 


4.2 Weak elastic coupling: Mott-lsing-to-elastic-Landau crossover 


For a relatively weak elastic coupling, on the other hand, there are two crossovers, see 
Fig.[^b), upon approaching the critical end point. First, a crossover at the Ginzburg 
scale occurs into a regime where the critical behavior is governed by non-trivial Ising 
exponents (blue shaded regime) until a second crossover occurs into the elastic critical 
regime (yellow shaded regime) described by Landau theory. 

In this section, we concentrate on the latter, Mott-Ising to elastic Landau crossover 
from the blue to yellow s had ed regime in Fig. [^b). For its description, we can adopt 
the approach of section (3.1), neglect the phonons and integrate out the Mott-Ising 
degrees of freedom. This yields an effective potential for the macroscopic strain singlet 


V{E) — + /sing(^j h + liE) 


(35) 


where we included already the corresponding singlet stress CTq. The function /sing 
is here governed by the Ising fixed-point and not known in closed form for general 
spatial dimension d. 

An extended Mott-Ising regime seems to be realized close to the Mott end point of 
the quasi two-dimensional charge-transfer salt k-(BEDT-TTF) 2 A 1 . Extensive investi- 
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gallons have demonstrated the critical properties in the vicinity of its end point is con¬ 
trolled by non-trivial critical behavior probably of Ising type |57l58l59l60l61l62l63l64j . 
It is suggestive that the reduced spatial dimensionality of the fluctuation in this ma¬ 
terial promotes the Ginzburg crossover and thus the Mott-Ising regime before the 
elastic crossover sets in. 

The available data of thermal expansion [63164) on k-(BEDT-TTF) 2 X allowed 
us in Ref. [53] to determine the parameters of the potential ( [3^ and to theoretically 
describe the expected second crossover to elastic criticality in this material. Motivated 
by its quasi-two-dimensional structure, we followed Bartosch et al. [Mj and used in 
Eq. ( j^ the free energy, /sing, of the two-dimensional Ising model, d = 2, in the 
asymptotic scaling regime, that can be computed numerically using the results of 
Ref. |7T| . Minimizing (35) with respect to E then yields the free energy E. 

Results from Ref. |33j for the T-dependence of the second derivative —d^E at the 
critical pressure Pc, which will govern the specific heat, are shown in Fig. I^a). Without 
the elastic coupling, 71 = 0, this derivative exhibits the logarithmic singularity of the 
2d Ising model (dotted line), that is however cutoff for finite 71 > 0 (solid line) and 
replaced by a pure mean-field jump [151 . The expected pressure dependence of the 
strain E for various temperatures is shown in Fig. [^b). Whereas for T < E the 
change of strain 6E as a function oi p — pc exhibits a first-order jump, it only exhibits 
a crossover for T > E. For T = E, finally, it behaves as SE \p — Pc\^^^ indicating 
a breakdown of Hooke’s law. 

Our quantitative analysis allowed us, in particular, to estimate the width of the 
elastic critical regime in k-(BEDT-TTF) 2 X along the pressure and temperature axis, 
Ap* and AE, respectively. 


Ap* = 45 bar, AT* = 2.5 K. (36) 

These values suggest that it should be experimentally feasible to detect the crossover 
from the Mott-Ising to the elastic critical regime, and thus investigate the change in 
universality class upon approaching the solid-solid end point of the Mott transition. 


5 Conclusions 

The characteristic feature of crystals is their shear rigidity which distinguishes solids in 
particular from liquids or gases. This shear rigidity and the corresponding long-range 
shear forces can in certain cases fundamentally alter the critical behavior close to a 
second-order phase transition. An especially drastic example is the elastic instability 
associated with a divergence of the bulk compressibility close to a solid-solid critical 
end point discussed in section [2?2| While the bulk compressibility diverges the phonons 
remain non-critical because the shear stiffness ensures that the phonon velocities all 
remain finite. As a result, the critical behavior of solid-solid end points do not possess 
a diverging correlation length and are described by true mean-field behavior. This is 
to be contrasted with the properties of the liquid-gas end point where the diverging 
compressibility implies the vanishing of the sound velocity, that in turn leads to a 
diverging correlation length and to non-trivial critical behavior belonging to the Ising 
universality class. 

In this paper we have discussed cases where the shear rigidity of solids strongly 
affects quantum phase transitions. Besides genuine elastic transitions, this physics is of 
importance for phase transitions that are primarily driven by magnetic, ferroelectric 
or electronic degrees of freedom. A first such case are critical end points of ’’liquid- 
gas type” described by an Ising order parameter. In such a case, the coupling of the 
order parameter to the strain field of the crystal will change the universality class 
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and critical elasticity prevails close to the critical point. We discussed two examples, 
the metamagnetic quantum critical end point in section and the Mott end point 
at finite temperature in section both of which are governed by the properties of 
solid-solid end points. 

Elastic quantum phase transitions associated with a change of crystal symme¬ 
try were discussed in section |2.1| In case of spontaneous crystal symmetry-breaking, 
phonon velocities in certain crystallographic directions vanish resulting in critical 
thermodynamics that is however still governed by a Gaussian fixed-point. We dis¬ 
cussed the corresponding quantum critical signatures and the resulting divergence of 
the Griineisen parameter. This quantum critical elasticity might be relevant, for ex¬ 
ample, for the pnictide superconductors at its continuous tetragonal-to-orthorhombic 
transition close to optimal doping [32]. 

Studies of the thermal expansion, compressibility and of ultrasound are especially 
useful to address the physics discussed in this review. With the help of such measure¬ 
ments, the strength of the elastic coupling can be determined and the crossover to the 
regime governed by critical elasticity can be quantitatively estimated, for a specific 
example see section [4. 2 1 Furthermore, they allow to identify directly a breakdown of 
Hooke’s law or the vanishing of sound velocities that are characteristic features of 
elastic phase transition. 


This is a review on work (Refs. | 32I33I34I41I48 |1 that was supported by the DFG through 
FOR 960. We acknowledge a collaboration with I. Paul and L. Bartosch on topics presented 
here. 
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